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A theory of sea scatter at large 

INCIDENT ANGLES 


A« K, Fung 
H# L, Chan 


ABSTRACT 


A theory for sea scatter at large incident angles (6^30®) is developed using 
q two'scale roughness model. The assumed small scale waves are to satisfy the small 
perturbation assumptions and the large scale waves to satisfy the physical optics 
approximations. Measured sea surface slope density and sea spectra reported by 
oceanographers are incorporated Into the theory to explain effects of incident angle, . 
polarization, frequency, wind speed and anisotropic characteristics of the sea surface. 
It is observed that the increase of the backscattering coefficients with the wind 
is due primarily to the growth of the seo spectrum and, to a lesser extent, to the inter 
action between the tvvo scales of reuyldicss. This iiMeiaciion effect is also the cause 
of the shift of the minimum of the scattering coefficient around the crosswind direction 
towards the downwind direction. The difference between the upwind and the crosswind 
observations Is the result of the anisotropic characteristics in the sea spectrum, while 
the difference between the upwind and dpwnwind observations is the consequence of 
the skewness in the slope probability density function of the large scale waves. Com- 
parison with some experimental data shows satisfactory agreements. 
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1,0 INTRODUCTION 


It has been generally agreed^ ^ thaf the scaftering properties of the sea 
surface can be explained reasonably well in terms of a two-scalo roughness model 
where the small scale waves are assumed to satisfy the small perturbation assumptions 
while the largo scale waves are to satisfy the Kirchhoff approximation. The basic 
approach is to use the first order result of the small perturbation method to nJompute 
the scattering coefficient due to small scale waves and to account for the tilting 
effect of the large scale waves by averaging this scattering coefficient over the slope 
distribution of the large scale waves. This averaged scattering coefficient is such 
that for near grazing incidence the vertically polarized scattering coefficient so 
computed is always larger than the corresponding horizontally polarized coefficient. 
Since recent experimental data^ indicate that there are cases where the converse is 
true, it is clear that such a theory may not be opplicable near grazing incidence. 
Hov/ovor, for incident angles between 30° and 80° the theory is expected to be valid, 

2.0 THE SCATTER TH60RY 

1 0 A 

Many authors have shown that for a slightly rough surface which satisfies 
the small perturbation assumption, the backscattering coefficient is given by (Figure 1) 

(0/0) = 2 K^°‘]^lcCpp|^W(0,0) (1) 

« 

where for horizontal polarization, p=h and 

Rj^ is the Fresnel reflection coefficient for horizontal polarization and for vertical 
polarization, p=v and 


where Rv and Tv are fhe Fresnel reflocHon and transmission coefficients for vertical 
polorization. In a , k Is the wave number in air; k‘ is the wave number in sea water 
and 0 is the angle oiT incidence. W (0,0) is the nomialized anisotropic sea spectrum 
and is the variance of the small scale sea waves* To include the tilting effect of 
the large scale waves, it Is necessary 'to average over the slope distribution 
(as viewed by the receiver) of the large scale waves. Thus, the incident angle In (1) 
becomes the local incident angle (denoted by 0' in (2)) and the averaged scattering 
coefficient is O^igure 1) 


05 03 


"■pp (Q,'^)=y f o-pp(0',0) P(, (&’, Zy')dZxdZy 


( 2 ) 


-05 -cotO . 

Note that the local Incident angle 0' is a function of the incident angle 0 and the 
surface slopes, Zk and Zy, In (2) Pq (Zx‘, Zy') is the slope distribution of the large 
scale waves as viewed at an incident angle 0 and is defined in the prime coordinates 
whose X* axis is parallel to the wind direction. It is assumed that the plane of in~ 
oidence is the XZ plane and that the angle between the x ~ axis and the x' axis is 
0 so that an upwind observation occurs when 0- 0, 

It is important to note that the probability of occurrence ot a slope on a given 
surface varies with the direction of observation. To illustrate this point consider a 
piece of one-dimensional surface depicted in Figure 2, The probability of occurrence 
For slope Zx 2 when viewed from an angle 0 is 


^ 2" 2 ° 


_ ^ 2"2 ° 


,Pe(Zxi.j= (fft i_e 


where a = zcos0 - xsin0 and n. are unit normal vectors to the surface. On the other 
hand, it becomes 


P (Zxz) == 




2"2 




^ 2"2 ^ 


L 


when viewed from ihe verHcal* Hence, 



I 

[ 2 ^ 2 ) = Iq (cos 0 + Zxp sin ^? ) P (Zx^) 


Since the total probability must bo one, the ratio 1/Lq can be determined by this 
condition. Similarly, for two-dimensional surface, the slope probability density function 
viewed from an angle 0 can be shown to relate to the slope probability density func- 
tion viewed from the vertical as 

P^ (Zx‘, Zy‘)= (I + Zx‘ tan0) P (Zx’, Zy') ( 3 ) 

Attempts to measure P(Zx*,Zy’) for the ocean surface were made by Cox and 
Munlc^ and for a slick sea surface (one which excludes small scale waves) it was 
reported that 


P{Zx',Zy‘) = 


F(Zx*,ZyO eXp ^ Zx'^ ^ Zy'^ 
2rrOy Og [ 2a^ ^ cTc \ 




* 


Each const•an^ in F{Zx',Zy‘) was reported by Cox and Munk to Ho within a 
range of values. Further Study by other ocRanographoR may lead to more refined 
estimates* For the purpose of illustration only one value Is chosen for each constant 
as follows; 

oro^=0.00740.78X10"^U 

ffc^-0*005+0.84X10“^U 

Cos=**0.05 C^o^O.Sd 
Co4c:0,26 

2 2 

In the expressions for (Ty and , U is the wind speed in motor per second at an 
altitude of 41 feet above the sea horizon. The relations between the surface slopes 
in the primed and unprlmed coordinates are as follows (see Figure 1); 

2x‘-2xcos0 +Zysin^^ 

Zy*= Zy cos0 « Zx sinS^ 


3.0 THE SEA SPECTRUM 

The spectrum for the small scale waves to be used in this paper is a modification 
of Pierson and Stacy's sea spectrum . The major change proposed by Pierson is that 
their spectrum expression for the capillary region has been replaced by the one reported 
by Mitsuyasu and Hondal®. The important characteristics to be noted in this spectral 
model are (1) the spectrum grows with the wind; (2) in the capillary region the larger 
the K number the faster Is the growth; and (3) this model is valid to about 38 knots . 

In accordance with Mitsuyasu and Honda, their model can be valid to friction velocity 
as high as 33m/sec at iO meters above the sea horizon. Thus, the theory may be valid' 
to higher wind speed than the 38 knots when operating frequency is such that only thejr 
portion of the sea spectrum is the significant contributor. Also, it should be noted that 
the definition of W(0,0 ) in (4) is not the same as the one used by Pierson and Stacy 


4 




In ihal W(6,0) is defined for all 0 between 0 and Sff# If Pierson and Stacy's sea 
spectrum is extended to cover the same angular range, it can bo shown that the com* 
ploto expression of the modified directional sea spectrum is os follows; 



5(K)0+a^cos2^) / K {'i) 

where higher order terms in 0 has been ignored; a^ is the only unknown parameter not 
yet specified; and where 


S(K)=S.(K), Kj,,<K<K. 


c /i^s_ a .... _r 0.74q^ 1 

[?U^J ' 

t 

S^CK) a K“^/- , K|< K<K2 == 0.359 

Sg (K) ~ 5^(Kg) (K/iCg/® , 1^2*^ ^ ^ * ^^2 

S4(K) = 0.87S(2it)'’’ ’ s i |fll<^Zl3.!?^ ,,, K3<K<Kj 

^ (gK+gKVl3.1769)<'^l*'’/2) 3 4 

Sj(K) = 1 .473 X K"^, K4<K< » 
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_L 


can bo found numortcally by sotHngs (K^) equal to (K^). 
U* « friction velocity, 

K^ = ( 13 , 1769)'/2 

Pj « 5,0 *^Iog^QU* 

s 0 ,684A + 4.28 X 1 0"^ v/ -4 .43 x 1 0**^ 

U(U*) = (U*/0,4) In (z/zj cm/soc 

a = 4,05 X 10 ^ 
g “ 980 cm/soc^ 

^*m“ ^2 cm/sec 


4,0 THEORETICAL RESULTS AND COMPARISON WITH EXPERIMENTAL DATA 


Upon combining (l) (3) and (4) and substituting into (2), we obtain for the 
bachscottered case 



S{2k sin 5V) (1 + Zx‘ tan 0 ) p (Zx‘, Zy') dZxdZy 



In (5) the unknown constant a^ can be estimated by plotting versus p curves and 
comparing these curves with experimental data. In all Figures 3-8 wind speeds are 
computed at 19,5m above the sea horizon. In Figure 3 and 4 ct° versus d curves are 
plotted for two different wind speeds at 6 * 30° and 6 = 60° and for various choices 
of a^. At 0 = 30° only vertically polarized curves are shown since the difference in 
a° due to polarization is small , At 60° the difference in cr° becomes significant and 
hence u for both polarizations are shown. The data points in Figure 3 and 4 were 
obtained from reference 11 . 
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In Figure 4a, the dola Indfcalcd by circles v/ero taken for 0 between 0 ® and 
180® and the squares denote data for d between 180® and 840®. The squaro-dqta 
were folded back under symmetry assumptions to compare with the circle data. It 
follows from 3o and 4a that a^-0,3 at 12 knots and a^J:;0.40 at 30 knots. It Is also 
noted from those curves that two local maxima occur around the upwind and the 
downwind directions and a local minimum occurs near the crosswind direction* In 
addition, the location of cr° minimum tends to shift towards the downwind direction 
as v/Ind speed or Incident angle Increases, This effect is most significant for horizontal 
polarization. The fact that or® upwind Is larger than a® downwind is due to the skew- 
ness in the slope density function, P 0 (Zx’,Zy‘). 

In Figure 4b and 4c 0 ® versus 0 curves using only the first order small per 
turbation theory are also plotted. Results correspond to the case when there is no 
tilting effect due to large scale waves. Comparison between these curves and those 
. obtained by using (5) Indicates that the use only of the first order small perturbation 
theory fails to account for (I) thn difference between the upwind and the downwind 
scattering coefficients (ii) the possible shift of the minimum of 0 ® towards downwind 
and (iii) the true 0 ® values around the crosswind direction which are higher duo to 
the interaction between the two scales of roughness. 

With a^ chosen to be 0.3, 0,4 and 0,50 at wind speeds at 6 , 2 , 12,6 and 18 
m/sec respectively, 0 ® versus wind speed curves are computed for two different 
frequencies, polarizations and incident angles as shown in Figure 5 through Figure 8 , 

In general, upwind scattering coefficient, 0 y, is found to increase foster with wind 
speed than the downv/ind coefficient, 0 p , which in turn increases faster than the cross- 
wind coefficient, 0 *^, All three scattering coefficients are found to increase faster 
with the wind speed as either the incident frequency or the incident angle increases* 

The rate of increase of the scattering coefficients with wind speed is larger for hori- 
zontal polarization than for vertical polarization. For incident angles beyond 60® 
the level differences between these coefficients also Increases with wind speed* When 
the theoretical 0 ® versus wind speed result is fitted over the wind speed range from 12,5 
knots to 25 ,4 knots by the equation 0 ® = AU where A is a constant; U is wind speed and 
y is the wind speed power coefficient, a table for V can be constructed for different fre- 
quencies, directions, polarizations, and incident angles. Such a table for two different 
frequencies Is given below, 

7 


§ 




theoretical v) 



8,9 GHz 

13.9 GHz 

DIRECTION 

POLARIZATION 

INCIDENT ANGLES 

INCIDENT ANGLES 



o 

o 

50° 

60° 

40° 

50° 

60° 

U 

W 

1,40 

1.48 

1.55 

1.63 

1*80 

1.81 

D 


1.38 

1.45 

1,49 

1.6 

1.75 

1.75 

C 


1.07 

1.16 

1.2 

1.3 

1.3 

1,45 

u 

HH 

am 


Hi 

1.70 

1.93 

2.01 

D 


Ib9 


Ira 

1.65 

1,82 

1.83 

C 


IBi 


H 

1.34 

1,43 

1.64 


To sop fha^ the wind depohdonco predicted from the present theory is in 
general ogrcomont with experimental data, a table showing the V values reported 
by Jones, ot end Moore, R. K., et al^^ is reproduced below. Comparison 
between these tables shows that general agreement ts, indeed, obtained for the 13,9 
GHz data, 



POIARJZAIION 

: y- 

mua.9citTi 

AAFS(I3.9GM/J 

SKYIAS S’I?? ()S.9CHit 

incicent ancle 

tNCI3.N1 ANGLE 

INCIOENT ANGLE 


so® 

AC.'* 

i0» 

AO® 


u 

W 



V.77 

1.66 

1.63 

1.62 

0 * 


o.co 



1,55 



c 



HI 

mm 

1.51 



mM 

HK 

0.0/ 

■SI 

t.w 

1.W 

1.63 

S.ll 



1,«. 


1.97 

1.96 



Bi 




H 


t.xs 
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5.0 CONCLUSION 


A grov/lng soa spoctrum^'^^ dnd a two-scolo roughnoss modol^"^ recommended 
by many Investigators for studying rqHor sea scatter have been extended to Include 
anisotropic characteristics of the sea surface. From such a theory, the following 
observations may bo made: 

(1) In the c® versus 0 ‘curves the local minimum occurs around the crosswind 
difoctlon and tends to shift towards the downwind side as the Incident angle 
Increoses* This is particularly true for horizontal polarization, 

(2) Stronger wind dependence occurs at larger incident angles end higher 
frequencies, 

(3) The wind depondonco of c® upwind Is stronger than cr® downwind which 
in turn is stronger than O'® crosswind # 

(4) Major factors contributing to anisotropy are the anisotropic sea spectrum 
and the nonzero mean slope In the large-scale- wave slope density function 
defined with respect to the plane perpendicular to the look direction. 
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ASPECT ANGLE (DEGREES) 

Figure 4b Azimuthal Dependence of o° at 6 =60° ,11=25^0 Knots 
and 13.9 GHz for Various Choices of for KH. 
The first order small perturbation results are 
also shovm tor comparison. 




ASPECT AMGLE (DEGREES) ■ 

Figure 4c Azimuthal Dapendence of a° at 0 =60® ,U=25.0 Knots 
and 13.9 GHz for Various Choices of 3o forW, 
The first order small perturbation results ai^e 
also shown for cemparison. ‘ 
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